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TECHNIGAL NOTE 3569 


CQMERESSIBLE LAMINAR BOUWEARY LAYER AND HEAT TRANSEER 
FOR UNSTEADY MOTIONS OF A FLAT PLATE 
By Simon Ostrada 


SUMMARY 

The laminar con^jressihle hotindary layer and heat transfer over an 
isothermal semi-infinite flat plate moving with a time -dependent veloc- 
ity has "been analyzed. First-order deviations frcan the quasi-steady ve- 
locity and ten^jerature profiles and boundary- layer characteristics have 
been cooxputed. For a plate oscillating about a steady velocity, it is 
sho™ that the Tnaxima. of skin- friction coefficient and local heat- 
transfer rate are out of phase with the plate velocity j the skin fric- 
tion leads by angles not exceeding 45° for peamissible values of the 
frequency parameter, idiereas the heat transfer is almost in’ phase with 
the plate velocity for very smal 1 Mach numbers but depends significantly 
on the Mach mmiber, plate to stream temperature ratio, and frequency for 
higher- speed flows. 


INTRODUCTION 

Until recently, studies of tinsteady laminar boundary layer were 
limited either to the early s-bages of the motion (i.e,, to the tran- 
sient state) or to oscillatory motions without a mean flow. The fluid, 
furthermore, 'vra.s assumed to be inccanpressible. More detailed investi- 
gations were not made, because it ’was felt that ’the boundary-layer 
grovrbh occurred in so short a period of time that, for engineering pur- 
poses, the flow could be assumed steady. However, in many present-day 
applications, consideration must be given to the unsteady effects for 
long periods of time and for high-speed flows in which ccsopresslbility 
is important. For example, the skin friction and heat transfer of the 
usual rocket missile must be regarded as unsteady for its entire flight 
because the flight speed varies continuously over the entire trajectory. 
Other cases of importance in this regard include blades rotating in non- 
uniform air streams, unsteady nozzle flow, and oscillating wings. 
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Accordingly, the unsteady lami nar ccmpressihle boundary layer over 
an insulated, surface uas analyzed, in reference 1. The develojpment 
therein is for continuoiis time-dependent velocities of the hody, and. 
universal functions are presented, from which the deviations of the ve- 
locity and. temperature profiles from the quasi-steady stated can he 
determined. 

To determine the effects of free-stieam fluctuations on hoth skin 
friction and heat transfer, the case of two-dimensional inccnrpressihle 
flow about a fixed cylindrical body is treated in reference 2 by an in- 
tegral method. 

In the present paper, prepared at the KAGA Lewis laboratory, con- 
sideration is given to the laminar compressible boundary layer and heat 
transfer over a semi-infinite flat plate maintained at a uniform (both 
temporally and spatially) tenperature and moving with a continuous but 
othervd.se arbitrary time-dependent velocity. The solutions are obtained 
as a series about the quasi-steady state. The analysis presented here- 
in, therefore, extends the results of reference 1 to include the effects 
of heat transfer. The present study also represents a more exact treat- 
ment including corapressibility effects of a similar^ problem treated in 
reference 2. 


ANALYSIS 
Basic Equations 

Consideration is given herein to the Jaminar flow and heat transfer 
about an isothermal semi-infinite flat plate in rectilinear flight 
through still air; the flight speed is to be time-dependent. For this 
flow it is presumed that the Prandtl boundary-layer assumptions are 
valid, and, in particular, the pressure is constant throu^out the 
fluid. If it is further assumed that the Prandtl number and specific 


^l-Jben the boundary layer at einy instant is that appropriate to 
steady flow at the instantaneous value of the stream velocity, the 
flow is said to be quasi steady. 

^The class of bodies considered in reference 2 is more general 
than that considered herein. 
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heat are constant^ the eq.-uations governing the flow and heat transfer 
in a cc» 5 )ressihle viscous fluid are 



\diere the symbols are defined in appendix A. 

Q?hese equations are written for a rectangular coordinate system 
vdilch is stationary in the fluid. 13ie plate,, therefore, moves ’t/lth a 
velocity U(t) in the negative x-dlrection, and y is measured normal 
to the plate (see fig. 1(a)). The origin of coordinates is taken to 
be the leading-edge location at t = 0. 

These equations may be written in another coordinate system fixed 
with reference to the plate, with the origin at the leading edge (see 
fig. 1(b)). The appropriate transfoimations for this change of coordi- 
nates are 


u=u + U^ v=v 

rt 


X = X + 


r u .t; 

'^0 


y = y; t = t 


Equation (l) then becomes 


35 + S + I? ‘ ° 

/Su ^ Su'S dU , S / Su^ 






be Be 

\3t + “ 35 + 


1 b f Be\ p /Bu\ 
^ V ^ Cp \5yJ 


Er 

pe = const 




(2) 
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The -vaxiation of viscosity with ten 5 )eratTJre is approximated hy 

P„v„C -f- (3) 

CK> CD 

as is discussed in reference 3. The constant C is obtained by match- 
ing equation (3) with the Sutherland formula at some appropriate point. 
If the matching is done at the toJ.1 , 


G =* 



(4) 


The momentum equation can be made independent of the energy equation by 
means of the transformation 


Y = 


dyj X = x; T = t 


( 5 ) 


'0 

which is similar to that used in reference 4. Further simplification of 
the basic equations 'vriLU result from the introduction of a stream func- 
tion as in reference 1: 


_ Rjd ^ 3 ilr 

^ = T 3y ^ ^ 






> ( 6 ) 


p , \ Poo , Silr By , 




and the definition 


9 = 


e - ©CO 

9w- 000 


Application of equations (4) to (?) to equations (s) yields 

^YT VxY “ %^YY ^ *^oo^YYY 


( 7 ) 


(8a) 


Gv„ Cv^ 

®T ^Y®X " "*^X®Y W ®YY Cp(0^ - 0^) ^^YY^ 


(8b) 
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I!he appropriate boundary conditions on ijr are 

^ ijfy(0,Y.T) = U(T) (9a) 

t(X,0,T) = 0 (9b) 

For the case of an isothermal surface as considered herein, the boundary 
conditions on 0 are 

0(X,«o,T) = 0(O,Y,T) = 0 (lOa) 

0(X,O,T) =. 1 (lOb) 

It should be noted that Initial conditions “would have to be added to 
these boundary conditions if the early stages of the motion are to be 
described properly. Hovrever, in the present problem it is assumed 
that sufficient time has elapsed so that the initial conditions no 
longer affect the flow. 


Solutions 

The solutions of the boimdary- value problem (specified by eqs. (s) 
to (lO)) describing the flow and heat transfer .of an isothermal semi- 
infinite flat plate traveling in a compressible viscous medium with a 
speed that may vary with time in a differentiable but otherwise arbi- 
trary manner will be obtained by a method similar to that of reference 
1. That method is appropriately modified herein to include the effects 
of heat transfer. 

Parameters . - It is desirable to determine the governing parameters 
before attenpting to obtain a solution. Eeference 1 noted that the di- 
mensionless parameters 

XU’. X^U”. . x^u(^). 

~~ 5 ~’ — — > ’ ' ' 

IP U* U“+l 

can be formed from the coordinate along the surface and the stream ve- 
locity and its derivatives. The Reynolds and Mach numbers are, of 
course, also pertinent. Physically, the quantities (U) represent the 
ratio of the time required for a change of some physical guantlty (e.g., 
velocity) at the boundary- layer edge to diffuse throu^out the layer to 
the time that is characteristic of a variation of the free stream at the 
boundary-layer edge. Hence, the quantities (ll) are a measure of the 
pronptness with which the boundary layer responds to inpressed changes. 
If the quantities (ll) are very small, the flow can be considered to be 
quasi steady; that is, the boundary layer at any instant is that appro- 
priate to steady flow At the instantaneous -value of the stream -veloci-ty. 
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The significance of the parameters (U) can he made clearer^ per- 
haps, hy consideration of a special case- For uniform acceleration, 
for exan 5 >le, U = ATj and the quantities (ll) reduce to a single parame- 
ter ^ ^ x/AT^, -which is equivalent to the ratio of the distance afl; of 
the leading edge to the dis-tance -traveled hy the plate. For ^ » 1, 
the effect of the leading edge has not yet heen felt at the station X; 
and reference 1 points out that the solution for this condition corre- 
sponds to the initial motion. If ^ « 1, the growth of the ho-undary 
layer -wi-th X must he considered; and reference 1 sho-vrs that for this 
case the flo-w- is nearly quasi steady- Returning to the earlier inter- 
pre-tation of (ll), it is clear that, if the velocity is increasing, the 
houndary Isyer at a fixed point on the plate hecomes progressively thin- 
ner and, hence, responds more quickly to changes at its outer edges. 

This trend tends to establish quasi-steady flo-w. 

Simplification for nearly quasi- steady flow-s . - In accord \'rlth ref- 
erence 1, if consideration is restricted to a s-tage of motion where the 
flow is nearly quasi steady^, the stream function may he defined as 

Vc^jiX f(a,^Q,^^,^2, . . .) (12) 


Y I U 

where cr - aItt — ^ and are functions of X, T that characterize 

the departure from the quasi- s-teady state. The explicit definitions of 
the are determined in the coxirse of the analysis. The dimension- 

less tenperatxire difference cannot he expressed in the lasiial manner (see 
ref. 3 ) in teamis of two functions, one of -tdilch solves the thermcaneter 
(i.e., insulated plate) problem and the other corrects for the heat 
transfer, since the adiabatic -wall tenperature is a function of both X 
and T herein. Therefore, let 

9 5 h(<r,£j + ( 1 =) 


^The motion in the early s-tages, that is, before the flow hecomes 
quasi steady, is treated for several configurations in reference 5. For 
this case, the Inertia terms and the x-dependence vanish. 
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Substitution of equations (l2) and (13) into equations (s) yields 


^oaa ~ 






0*= 


it=0 / 


f + 2 ^ Of 
^ ^ y2 ^0-cr 


w w 

- 2X t„a E \ ?nX + ‘ i ^ 5:ffl 

n=0 n=0 


(14) 


iW) 


^ U ^ ■*■ ^ Hn ^nX “ ^ ^nX 


XpmO 




(15) 


s„<, + PrtBa + ^ - 4Er| 


I ^ 2 ^ ® i ^ ®£n if“ ^ 

_ \ z>*0 n»0 / 


(16) 


Equations (15) and (l6) are coirect only if the quantity Cv„ is con- 
stemt. For the case considered herein (constant 0^)^ equation (4) 
shows that Cv„ is indeed constant. !Phe boundary conditions (eqs. 

(9) and (lO)), in terms of the definitions (l2) and (13), are 


■n) - 2; 

fp(0,^^) f(0,^^) == 0 

(17) 

= 0; 

n(o,Cj == 1 

(18) 

s(»^^n) = 0; 

s(0,^J = 0 

(19) 


Since nearly quasi -steady flows are assumed here^ no Initial conditions 
need be specified. 


Equations (l4) to (16) can be made self-consistent (i.e., functions 
of a and only) by defining 

^ _ XU‘ V X^U” ^ X^U’" 

&0= ^1= ^2= • • • 

These definitions correspond to the quantities (ll); and_, since the ra- 
tio of diffusion time to the characteristic free-stream variation time 
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and, hence, the given hy (20) are assimied smal.l relative to unity, 

the functions f, h, and s may he e;j 5 )anded as follows: 


f + 






H(ct) + ^QbQ(cr) + + • • • + . . . + ^o^l^l^'^^ **■ 


S(cr) + ^qSq(ct) +^iSj^(ct) + . • . + . • - + 


( 21 ) 


( 22 ) 


(23) 


A disciission of the magnitudes of the is presented in reference 1. 

It is there pointed out that in practical situations the requirement of 
sman is commonly met; and, furthermore, the generally form a 

diminishing sequence, provided that U(t) is a differentiable function. 
The li mitations for small for various specific U(t) are also dis- 

cussed in reference 1. 

Siibstituting equations (2l) to (23) into equations (14) to (19) 
and collecting coefficients of the various powers and products of the 
yield in part (up to order ^^) 



F«» + FF" = 0 

(24a) 

f ttt . Jljf'I _ 
0 0 

2F'f^ + 3F”fQ = - 4(2 - F') + 2 oF" 

(24b) 

1 

+ Ff^ - 4F'f| + 5F"f^ = 4f^ 

(24c) 


F«(<*.) = 2; F*(0) = F(o) = 0 (25a) 

fi(») = f^(0) = f^(0) = 0 (25b) 

H" + ErFH’ =• 0 (26a) 

h" + PrFh^ - 2PrF'h^ = FrH* (2a - Sf^) (26b) 
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hj + ErKij' - 4PrF'h^ = Et(411q - 5H’f^) (26c) 

H(a.) = 0; H(0) = 1 (27a) 

i4,(«) = hj^(o) = 0 (27b) 

S" + ErFS' + ^ (F")2 0 (28a) 

s" + PrFs^ - 2 PtF'Sq = Er(2oS‘ + 83 - F"fJJ - 33 ’f^) (28b) 

+ ErFsj^ - 4PrF's^ = Pt(46q - F"f^ - 53'f3_) (28c) 

3(«) = 3(0) = 0 (29a) 

Sn(®) = Sq( 0) = 0 (29b) 


Equations (24a), (25a), (26a), (27a), (28a), and (29a) are the 
equations appropriate for steady flow, -where, of course, F is -the 
Blasius function, -where H describes the ten 5 )erattire dis-feribution 
for an isothermal plate if the dissipation is neglected, and -vdiere 

3(c) => R(ct) - K(0) H(a) (30) 

The R(cr) is the solution for an insulated plate, and R(0) is the re- 
cQ-very factor. The functions F, H, and R are -tabulated in numerous 
papers; for exanq)le, F and H can be found in reference 6, and R is 
given in reference 1. 3ince the cited eq-uations represent steady flow, 
equations (2l) to (23), accordingly, indicate -that for small the 

laminar compressible boundary layer on an iso-thermal flat plate is 
nearly quasi steady -v/ith respect to both -the -veloci-ty and tempera-ture 
distributions. Hence, the first-order deviations from -the quasi-steady 
case v/ill be determined by solving equations (24b), (24c), (26b), (26c), 
(28b), and (28c) v/i-th their associated boundary conditions. 

3olutions of momentum equations . - 3ince the momentum equation ’tra.s 
made independent of the energy equation by -the -transformation given by 
equation (5), -the solutions of equations (24) and (25) are not altered 
by a change in -the theamial boundary condi-tions. Therefore, the func- 
tions and their deri-vati-ves for the case of an iso-thermal 

plate considered herein are identical to those for an insulated plate 



as treated In reference 1. The function F and its derivatives^ as has been stated, are tabu- 
lated for Pr a 0.72 in reference 6, eind the functions fQ, fl, and their derivatives were de- 
termined for Pr = 0.72 by direct numerical integration in reference 1 and are presented 
therein. 


Solutions of energy equations . - The solution of equations (26a) and (27a) is known and, 
as has been stated, is tabulated in reference 6. The function H is presented in table l(a) 
and figure 2(a) of this r^ort for cooipleteness . The function S, presented in table l(b) and 
figure 2(b), is determined frcan equation (30) using the known H and R, the latter being tab- 
ulated for Pr 0.72 in reference 1. The ranalning energy equations and their associated 
boundary conditions (eqs. (26b), (26c), and (27b), and (28b), (28c), and (29b)) are solved for 
Pr » 0.72 by the nimerical Integration methc^ described in detail in appendix 0 of reference 1. 
An outline of the method is presented herein in appendix B. The functions hQ and hj_ are 

presented for Pr = 0.72 in table l(a) and figure 2(a), and Sq and for Pr » 0.72 are 

given in table l(b) and flgvn*e 2(b), 


UNSTEAOT BOUMDAHY-LAim CHARACTERISTICS 
Velocity amd Tengierature Profiles 
The velocity and temperature distributions can be obtained from 


H « i|p«(a) + ^Qf^(a) + + . . ^ 


0-00, U2 

^ a H + 

-• 


2Cp(0v - ej S + ^0 ^ + 2cp(0v - 0«) 


+ ^2. 


hn + 




•1 2Cp(0^ -■ 0-J 


The relation between y and cr is 


/CV.3 


( 7 +- 


H dcy+J 
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The functions F and fn are independent of the Mach numher. 
Hence^ the velocity profile is thickened (or thinned) because of com- 
pressihility. Tt^e specific effect can he determined from the preceding 
relation between y and a and depends on the temperature ratio and 
the Mach number. The unsteady effects are associated with the terms 
containing the • The functions hQ_, hj_^ Sq, and S]_ associ- 
ated with the temperature profile are presented in figure 2. 


Skin Friction 


The local skin- friction coefficient may be written as 


C4> = 


'•w 

^ P 
2 


1 

2 


& 

UX 


F"(0) + + ^ifi(o) + 



where = [^|-i(^/Sy)J-w is the wall shear stress. Substituting the 
values of F"(0), fQ(0), and fp(o) for Pr =3 0.72 given in references 
1 and 6 yields 

1 + 2.555 Co - ^1 + • • •] 

Equation (3l) has the same form as its comterpart for an insulated sur- 
face as given in reference 1, btrb the skin friction does differ for the 
case of heat transfer. The effect of heat transfer is accounted for by 
the constant C (see eq. (3)). From equation (4) it can be seen that C 
will be altered, because dy is the isothermal value specified herein 
rather than the adiabatic value used in reference 1. 


Cj = 0.6641 





The leading term on the right of equation (3l) is the quasi-steady 
value of the coefficient. As would be expected, positive acceleration 
leads to larger values of skin friction than the quasi- steady value. 


Displacement Thickness 


The displacement thickness 5* is defined as 


6 ^ 



_P_ 

p-u; 


dy 





or^ uae of equations (2)^ (s)^ (s), and (l2), 




CO 



Introducing equations (?) and (l3) gives 


.1^ F 


d 


12 + 


2(©w - ®-) 


h + (r - i)m^s - 


dcr 


0 


or^ using equations (2l) to (23); 




P)+ (y- l)H^ 






-CoM + (r-l)n‘^a(/la+ ^f^''l^^^ ^hgaa+ . . . +Ci -fi(») + (T-l)Nj,^ild« + ^i^2^^^^lde+ . 


Using the fact (ref. l) that 

CT— r 5 _ - 

i.m *ij^^ |i+ o.1673(y - i)«^+ i.iasa - 


11m (P - 2cj) *t - 1.721 and the Infojmiation given in table I; 

a-*co 

r r ,__Tl 

0.6919 qj|l+ 0.4M9 (y - 1)n|i - 0.0e«88 + 0.#Bti Cl|l+ 0-2 SBO(y - 1)hS + 0-5285 


(32) 


!The first three terms on the rl^t correspond to the ccii 5 )reB Bible quBsi-steady result; as can be 
verified for a given Mach number and tenroerature in reference 6. The effects of heat transfer 
not only affect the quasi-steady results through G and (0^ - but also are of impor- 

tance in the deviations frcm the quasi steady. For example; if 9-^ < 9„, a positive acceler- 
ation will lead to a thinner boundary layer than the quasi steady; as vas the case also for an 
insulated surface. However; if 9^ > 9^, the boundary layer will depend on the magnitude of the 

Mach number as well as the teinperature difference. 
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HociS-X RSr'bo of Hg Bf'tf lEr^snsf©!* 
The local rate of heat transfer is given by 


q. 3 -kl 


w 


Using the definition of the Prandtl number and equations ( 2 ), ( 3 ), (s), ( 7 )^ and (l3) in the 
above expression yields 

-r 


q w 




- 2 § “ 0 «) 
or, using equations (22) and (23), 






«i(0) + . . - 


Substituting the numerical values frcm table 1 yields for ih: = 0.72 

1-0.4106 


°p(«v - " 0.4240(r - X)J^ g-~ - 0.06923 £ojl+ 0.274s(r- ‘ °'*®“ ‘ 


CM) 


The compressible quasi-steady heat-transfer rate is given by the first two terms on tne 
right side of equation (33), as can be seen by ccmparing with equation (21) of refe;rence 6. 


The deviations ffcm the quasi-steady results depend again on the magnitudes of the temper- 
ature difference and Mach number. For example, for positive acceleration the heat-transfer rate 
from the plate is decreased if 0^ if imLght either be Increased or 

decreased. 
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Plate Oscillating About a Steady Velocity 


As an exan^le of tlie foregoing analysis, sTjppose that the plate os- 
cillates about a stea^ velocity as 


U =s Uq,( 1 + 6 sin (fit) 


(34) 


■where the amplitude of the velocity fluctuations e is to he mnal 1 rel- 
ative to unity. Substituting equation (34) into equations (3l) and (33) 
yields to order e 


and 


Cf^^ + 6C3^sin(mt + P]^) 


q =» <3 q|i + £C 2 Sin(cDt - (Pg) 


(35) 


(36) 


where c^^^ and qQ are the steady (corresponding to flight velocity 

Uq) local skin- friction coefficient and heat-transfer rate, respec- 
tively, and ■5'diere 


^ 1.5 + 3 


=. tan"^ 1.703^^j - 1 d- 


(37) 

(38) 


for p 4 

^ 2.120 


Cg - 


1 - 2. 1203 J 


2(1 - 0.4243) 


1 + 0.8560 - 2.25503 + 0.98833^ ^ ^ 

(1 - 2.1203)2 V^J L\i^/ J 


(39a) 


<Pg = tan' 


0.1385(1 + 0.27463)1 — 

-1< W 


1 - 2.1203 


1 + erl 


mi 


and, for B = , 

, ^ 2.120^ 


C? = 0.09775 — 


[‘•"[©‘I} 


(40a) 


(39h) 
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VThere 



(40b) 


(41) 


The quadrants of the angles (p]^ and <^2 can be determined from the 

respective signs of the ntimerator and denominator of the defining 
equations . 

It is significant to note from equations (34) to (36) that the max- 
ima of skin friction and heat transfer are not in phase with the maxima 
of the plate velocity; and^ furthermore, the amplitude of the fluctua- 
tions of these quantities depends on the frequency, as was also found in 
reference 2. 


Note that equations (35) and (36) give the deviations from the 
steady (or mean) conditions lattier than from the quasi steady as was 
the case in the discTission of the general analysis. The functions C^^ 

and C2 ('tdiich define the maximum deviations from the steady) and the 
phase angles and 93 are presented in figure 3. It can be seen 

from the definition of and figure 3(a) that for permissibly values 

of — the maxima in skin friction will always lead the plate velocity 

by angles not exceeding about 45°. The phase relations betireen the 
heat-transfer rate and the plate velocity depend essentially on the Mach 
manber, ratio of plate to stream ten^erature, and frequency. For Mach 
numbers near zero, ^3 beccsnes small, and the maxima of hyat transfer 
and plate velocity are nearly in phase in accord with the discussion in 
reference 2. 


From figure 3 it can be seen that the unsteady effects can be sig- 
nificant. It should be noted from equations (36) and (39) that 03 be- 
comes infinite >±ien the steady heat transfer approaches zero. The ac- 
tual heat-transfer rate given by equation (36) is, however, finite. The 

function C3 remains almost constant with — (see fig. 3(b)), so that 

the maxlmxmi anplltude of the heat-transfer fluctuations is essentially 
the quasi -steady value. 
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CONCLUDING EEMARKB 

The laminar con^jressihle botindary-layer flow and heat transfer over 
sin isothermal sem i . -infinite flat plate moving with a time-dependent ve- 
locity has been analyzed. The first-order deviations of the velocity 
eind ten 5 )erattire profiles, skin- friction coefficient, displacement thick- 
ness, and local heat-transfer rate from the quasi steady have been com- 
puted; the associated •universal functions are presented in tabular form 
for Pr = 0.72. Eelative to quasi-steady flow, positive acceleration 
results in larger skin friction, thinner boundary layers if the wall 
tenrpera-tnre is larger than the free-stream temperature, and either 
thicker or thinner boundary layers (depending on the Mach number ) if the 
sxirfa.ce is at a loxrer temperature than the stream. Positive accelera- 
tion results in lower heat-transfer rates from the plate if the surface 
temperature is greater than the free stream, ■whereas 'vrLth lower surface 
tempera'turas 'the heat-transfer rate -would be increased or decreased de- 
pending on the magnitudes of the temperature difference and the Mach 
number. Hence, the boundary-layer characteristics for an Isothermal 
surface can be considerably different from those for an insulated 
surface . 

Consideration of the particular case of a plate oscillating about 
a steady velocity sho^ied that the boundary-layer characteristics can be 
appreciably altered by the unsteady effects. Furthermore, tha skin 
friction and local heat -'transfer rates -were found to be out of phase 
■with -the plate velocity for peamissible -values of the frequency param- 
eter. The mawi.ma of skin friction lead the plate -valocity by amounts 
not greater than about 45°. For steady Mach nxmibers_near zero, the heat 
-transfer is almost in phase wi-th the plate velocity. At higher speeds 
the heat--transfer phase angle depends significantly on the steady Mach 
nxmiber and ratio of wall to free-s-tream tempera-ture in addition to the 
frequency. 


Le-wis Flight Prop-ulsion Laboratoiy 

National Ad-visory Committee for Aeronautics 
Cle-veland, Ohio, Axigust 10, 1955 
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APPENDIX A 
SYMBOLS 

A acceleration 

C constant defined by eq. (4) 

C]_^C 2 constants defined by eqs. (37) and (39), respectively 
Cj local skin- friction coefficient 

Cp specific beat at constant pressure 

P related to stream function for flat plate In steady flow 

f,f^ functions related to stream function for unsteady flat-plate 
flow, 1 =3 0, 1, 2, . . . 

H ten^jerature function related to steady flat-plate flow 

b,h^ functions related to temperature for unsteady flat-plate flow, 

1 = 0 , 1 , 2 , . . . 

k thermal conductivity coefficient- 

M Mach number 

Pr Prandtl mmiber 

q local heat-transfer rate 

R function related to tenqjerature profile for Insulated . flat plate 

In steady flow 

S function related to tenqoerature for steady flat-plate flow 

s,Sj^ functions related to ten^erature for unsteady flat-plate flow, 

1 = 0 , 1 , 2 , . . . 

T,t,t time 

U stream velocity In x- direction 

u, u velocity In X- and x-directlons, respectively 

v, v velocity In Y- and y-directlons , respectively 
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X,x,x coordinate along surface 

coordinate norma], to surface 
3 constant defined, "by eq., (il) 

T ratio of specific heats 

6* displacement thickness 

6 amplitude of velocity fluctuations 

dimensionless parameter, n = 0, 1, 2, . . . (eq. (20)) 

0 dimensionless ten^jerature difference 

0 ten^jerature 

p. absolute viscosity coefficient 

v„ kinematic viscosity coefficient outside boundary layer 

p density 

Y 

cr dimensionless coordinate, — 

<Pl,<p 2 phase angles defined by eqs. (38) and (40), respectively 
i[r stream function 

OQ frequency of velocity fluctuations 

Subscripts : 

w evaluation at wall (y = O) 

0 evaluation at a steady condition 

oo evaluation in stream (y -»■ <») 

Subscript notation for partial differentiation is used idaen con- 
venient. Primes denote ordinary differentiation. 
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APPENDIX B 


IKIEGEATION MEHHOD 


The second-order differential equations for the ten^Derature func- 
tions (eqs. (26b), (26c), (28b), and (28c)) and their associated bound- 
ary conditions (eqs. (27a), (27b), (29a), and (29b)) constitute a tvro- 
point boundary- value problem. The difficulty presented by the fact that 
the boundary conditions are given at two points (a = 0,®) is overcome by 
splitting the original two-point problem into two' single-point problems. 
To this end, each function Bq, h]_, Sq, and S]_ is written as the sum 

of two functions; for exanqple, for Bq, as 


ho= - 


1 1m 


4^^ 


(ct) 




o^oo h^'^''(cr) 




h^ Htj) + h^ ^(ct) 


(2). 


0 


,(D- 


where the (a) and h^^^(cr) are solutions of single-point problems. 


.( 1 ) 


The Bq •' satisfies the hcmogeneous equation iri.th lie specified Initial 


condition plvis an arbitrary initial condition 
the boundary condition at infinity, and the h, 


' ( 0 ) 


1 replacing 
satisfies the non- 


homogeneous equation and homogeneous initial conditions on the function 
and its first derivative. 


The scale of the variable a is then divided into equal intervals. 
Starting from the initial values of each part of each function - for ex- 

an^le, and - the values at successive points near 

0^0 can be obtained by expanding the unknown function in a Maclatorin 
series. Thus; the function and its derivatl'ves are knoim at a succes- 
sive number of points near a =« 0. A polynomial (for the highest deriv- 
ative) can then be matched to the known values and to the unknown value 
at the nex± point. The degree of the polynomial and size of the inter- 
val depend on the accuracy required. In this regard, it can be stated 
that the solutions of the energy eq^iations are more easily obtained than 
those of the momentum equations, since a second-degree polynomial is 
matched at three successive points for the temperature functions, where- 
as a fifth-degree polynomial was matched at six points in the solution 
of the momentimi equations (see ref. l). This simplification is consid- 
ered warranted, because the energy equations are of lower order than the 
momentum equations. However, because the curvatures of the sq and S]_ 
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functions were relatively greater tlian tlie 1 iq and h]_ functions (see 
fig. 2.), it was decided to check the accuracy of the three -point method 
described herein. Accordingly, sq and its derivatives \Texe obtained by 
five-point integration formulas, and those results for Sq differed by 

at most one in the fourth decimal place and by two in the fourth place 
for Sq from those presented herein. 

The polynanlal is then integrated to yield the successively lower- 
order derivatives at the unknown point in terms of the kno™ and imknown 
values of the highest derivative. The condition that the function and 
all its derivatives must satisfy the differential equation at the im- 
kno^m point serves to determine the highest-order derivative eind, hence, 
the function at that point. Thus, given the function at a successive 
number of points, the solution can be extended to the next point. This 
procediore is then continued over the entire range of cr. 
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TABLE I. - FUNCTIONS ASSOCIATED WITH TEMPERATURE PROFILE 
(a) Functions h. 


a 

H 

H' 





0 

1.0 

-0.5912 

0.00000 

0.04093 

0.00000 

0.25019 

.1 

.9409 

-.5911 

.00397 

.03737 

.02503 

.25039 

.2 

..8818 

-.5905 

.00732 

.02915 

.05004 

.24939 

.3 

.8228 

-.5887 

.00976 

-.01934 

.07400 

.24500 

.4 

.7641 

-.5853 

.01122 

.01014 

.09887 

.23539 

.5 

.7058 

-.5796 

.01186 

.00300 

.12167 

.21932 

.6 

.6483 

-.5713 

.01192 

-.00135 

.14251 

.19631 

.7 

.5917 

-.5599 

.01170 

-.00275 

.16070 

.16659 

.8 

.5364 

-.5452 

.01146 

-.00152 

.17562 

.13108 

.9 

.4827 

-.5269 

.01145 

.00169 

.18677 

-09127 

1.0 

.4311 

-.5050 

.01183 

.00603 

.19379 

.04908 

1.1 

.3819 

-.4797 

.01266 

.01057 

.19657 

.00660 

1.2 

.3353 

-.4512 

.01391 

.01448 

.19518 

-.03405 

1.3 

.2917 

-.4199 

.01551 

.01708 

.18989 

-.07092 

1.4 

.2513 

-.3865 

.01727 

.01796 

.18118 

-.10244 

1.5 

.2144 

-.3516 

.01904 

.01697 

.16963 

-.12749 

1.6 

.1811 

-.3161 

.02061 

.01423 

.15591 

- .14548 

1.7 

.1512 

-.2805 

.02184 

.01002 

.14076 

-.15633 

1.8 

.1249 

-.2458 

.02258 

.00479 

.12487 

-.16042 

1.9 

.1020 

-.2125 

.02278 

-.00082 

.10887 

-.15850 

2.0 

.0824 

-.1812 

,02242 

-.00631 

.09333 

-.15159 

2.1 

.0*657 

-.1524 

.02153 

-.01122 

.07868 

-.14086 

2.2 

.0518 

-.1264 

.02020 

-.01524 

.06523 

-.12749 

2.3 

.0403 

-.1034 

.01851 

-.01816 

.05321 

-.11264 

2.4 

.0310 

-.0834 

.01660 

-.01992 

.04272 

-.09730 

2.5 

.0236 

-.0663 

.01457 

-.02057 

.03374 

-.08227 

2.6 

.0176 

-.0520 

.01252 

-.02025 

.02623 

-.06817 

2.7 

.0131 

-.0401 

.01055 

-.01919 

.02006 

-.05538 

2.8 

.0095 

-.0306 

.00870 

>- 01760 

.01510 

-.04415 

2.9 

.0069 

-.0230 

.00703 

-.01564 

.01117 

-.03461 

3.0 

.0049 

-.0170 

.00557 

-.01351 

.00812 

-.02665 

3.1 

.0035 

-.0124 

.00432 

-.01136 

.00579 

-.02019 

3.2 

.0024 

-.0089 

.00329 

-.00933 

.00404 

-.01507 

3.3 

.0016 

-.0063 

.00246 

-.00748 

.00274 

-.01107 

3.4 

.0011 

-.0044 

.00179 

-.00588 

.00180 

-.00805 

3.5 

.0007 

-.0030 

.00127 

-.00452 

.00110 

-.00581 

3.6 

.0005 

-.0021 

.00087 

-.00342 

.00061 

-.00418 

3.7 

.0003 

-.0013 

.00057 

-.00254 

.00026 

-.00303 

3.8 

.0002 

-.0009 

.00036 

-.00184 

.00000 

-.00221 

3.9 

.0001 

-.0006 

.00020 

-.00134 

-.00019 

-.00168 

4.0 

.0001 

-.0004 

.00009 

-.00096 

-.00035 

-.00136 

4.1 

.0000 

-.0002 

.00001 

-.00068 

-.00048 

-.00116 



= 0.9692 

So 

= 0.04425 

dff = 0.3158 
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TABLE I. - Concluded. FUNCTIONS ASSOCIATED WITH TEMPERATURE PROFILE 

("b) Functions s. 


a 

5 

S ' 

"0 

-0 

"1 


0 

0.00000 

0.50134 

0.00000 

0.02248 

0.00000 

- 0.27540 

.1 

.04694 

.43775 

-.01217 

-.24852 

-.01863 

-.09863 

.2 

.08754 

.37395 

-.04648 

-.42273 

-.02027 

.06196 

.3 

.12176 

.30988 

-.09399 

-.51508 

- .00719 

.19431 

.4 

.14953 

.24593 

-.14724 

-.53969 

.01741 

.29163 

.5 

.17099 

.18287 

-.20013 

-.51025 

.04987 

.35127 

.6 

.18614 

.12174 

-.24794 

-.44003 

.08643 

.37400 

.7 

.19540 

.06368 

-.28721 

-.34181 

.12356 

.36311 

.8 

.19905 

.01006 

-.31577 

-.22761 

.15812 

.32417 

.9 

.19766 

-.03789 

- .33257 

-.10845 

.18770 

.26433 

1.0 

.19174 

-.07898 

-.33761 

.00616 

.21053 

.19060 

1.1 

.18205 

-.11227 

-.33176 

.10850 

.22561 

.11051 

1.2 

.16953 

-.13732 

-.31651 

.19305 

.23264 

.03094 

1.3 

.15491 

-.15411 

- .29384 

.25662 

.23200 

-.04225 

1.4 

.13905 

-.16293 

-.26591 

.29822 

.22455 

-.10467 

1.5 

.12261 

-.16462 

-.23490 

.31885 

.21152 

-.15350 

1.6 

.10625 

-.16019 

-.20276 

.32099 

.19434 

-.18760 

1.7 

.09070 

-.15159 

-.17119 

.30819 

.17447 

-.20734 

1.8 

.07615 

-.13899 

-.14149 

.28445 

.15329 

-.21424 

1.9 

.06297 

-.12436 

-.11454 

.25375 

.13199 

-.21036 

2.0 

.05125 

-.10883 

-.09085 

.21959 

.11150 

-.19835 

2.1 

.04119 

-.09328 

-.07065 

.18486 

.09250 

-.18082 

2.2 

.03259 

-.07843 

-.05384 

.15171 

.07544 

-.16012 

2.3 

.02548 

-.06475 

-.04021 

.12155 

.06051 

-.13820 

2.4 

.01961 

-.05254 

-.02941 

.09517 

.04777 

-.11657 

2.5 

.01485 

-.04196 

-.02104 

.07283 

.03714 

-.09628 

2.6 

.01121 

-.03295 

-.01471 

.05449 

.02844 

-.07797 

2.7 

.00822 

-.02558 

-.01003 

.03985 

.02146 

-.06202 

2.8 

.00605 

-.01945 

-.00665 

.02842 

. 0159.7 

-.04848 

2.9 

.00431 

-.01460 

-.00427 

.01979 

.01169 

-.03723 

3.0 

.00305 

-.01082 

-.00263 

.01339 

.00844 

-.02814 

3.1 

.00207 

-.00790 

- .00152 

.00878 

.00601 

-.02093 

3.2 

.00143 

-.00570 

-.00082 

.00555 

.00421 

-.01533 

3-3 

.00098 

-.00405 

-.00039 

.00336 

.00289 

-.01107 

3.4 

.00063 

-.00283 

-.00014 

.00193 

.00196 

-.00787 

3.5 

.00042 

-.00198 

.00002 

.00101 

.00130 

- 1.00551 

3.6 

.00021 

-.00129 

.00008 

.00046 

.00085 

-.00381 

3.7 

.00013 

-.00096 

.00011 

.00013 

.00053 

-.00259 

3.8 

.00004 

-.00060 

.00012 

-.00004 

. 0003f 

-.00176 

3.9 

.00002 

-.00037 

.00010 

-.00012 

.00017 

-.00117 

4.0 

.00005 

-.00023 

.00009 

-.00015 

.00008 

-.00077 

4.1 

0 

-.00019 

.00008 

-.00015 

.00002 

-.00051 


Jo' ^ " 

= 0-2879 

v/o 

= - 0,4440 

\ 


= 0.3095 
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(b) Functions s. 

Figure 2. - Functions associated with temperature profile 
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